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A relativistic  computer  simulation  demonstrates  that  the  application  of  a DC 
electric  field  enhances  the  excitation  of  plasma  waves  through  runaway  electron 
The  presence  of  an  elongated  tail  in  the  distribution  function  gives  rise  to 
acoustic-type  modes  which  interact  strongly  with  the  runaway  population.  The 
waves  produce  a strong  drag  which  reduces  the  runaway  momentum  Increase  from 
2/3  to  1/2  of  the  free  fall  value.  A threshold  DC  electric  field  Is  found 
above  which  there  is  no  collective  slowdown. 
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ABSTRACT 

— 

A relativistic  computer  simulation  demonstrates  that  the  application 
of  a DC  electric  field  enhances  the  excitation  of  plasma  waves  through 
runaway  electrons.  The  presence  of  an  elongated  tail  in  the  distribution 
function  gives  rise  to  acoustic-type  modes  which  interact  strongly  with 
the  runaway  population.  The  waves  produce  a strong  drag  which  reduces 
the  runaway  momentum  increase  from  2/3  to  1/2  of  the  free  fall  value.  A 
threshold  DC  electric  field  is  found  above  which  there  is  no  collective 
slowdown. 
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The  evolution  of  runaway  electrons  taking  into  consideration  the  self- 
consistent  excitation  of  large  amplitude  collective  oscillations  is  of 
considerable  interest1’2'3  at  the  present  time.  This  interest  has  been 
stimulated  by  recent  experimental  observations*4’ 5,6  in  tokamak  devices  of 
runaway  electrons  accompanied  by  the  spontaneous  emission  of  intense  radia- 
tion in  various*4’5  frequency  intervals.  The  study  of  runaway  electrons 
provides  a suitable  arena  to  test  and  develop  various  theories  of  collective 
interactions  in  plasmas  whose  applicability  can  extend  beyond  just  the  runaway 
problem.  Recent  theoretical  studies1  rely  on  assumptions  concerning  a flat 
momentum  distribution  parallel  to  the  applied  electric  field.  However, 
many  details  associated  with  the  wave-particle  interactions  that  generate 
and  alter  such  distributions  still  remain  to  be  illustrated.  One  approach 
previously  used  to  examine  some  of  the  relevant  physics  has  been  the  Fokker- 
Planck  type  of  simulation7’8  in  which  the  wave-particle  interactions  are 
assumed  to  be  described  by  quasilinear  theory.  The  latter  assumption  is 
also  characteristic  of  most  of  the  theoretical  treatments  of  the  problem. 

In  yet  another  work,  non-relativistic  particle  simulations  without  DC  field 
and  starting  from  a flat  distribution  have  been  performed8  to  verify  that 
collective  effects  could  lead  to  a positive  slope. 

This  study  examines  some  of  the  wave-particle  interactions  and  the 
nature  of  the  collective  modes  present  in  the  evolution  of  runaway  electrons, 
employing  a full  particle  computer  simulation.  The  major  effects  found  are: 
(1)  the  runaway  electrons  exhibit  large  excursions  in  momentum  space  not 
amenable  to  a quasilinear  discription,  (2)  the  momentum  distribution  function 
develops  an  elongated  tail  which  gives  rise  to  modes  with  frequencies  lower 
than  the  plasma  frequency,9  (3)  the  excitation  of  large  amplitude  waves 
reduces  the  rate  of  increase  of  the  average  runaway  momentum  from  1/2  to  2/3 
of  the  free  fall  value. 
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A one  dimensional  relativistic  electrostatic  particle  code  is  used. 

The  system  is  initialized  with  a low  density  cold  electron  beam  and  a 

Maxwellian  background  distribution.  The  external  X electric  field  is 

applied  selectively  to  fast  particles  with  momentum  above  a set  value, 

chosen  smaller  than  the  initial  beam  momentum  and  larger  than  the  initial 

average  momentum  associated  with  the  thermal  velocity  of  the  background 

particles  vTQ.  This  technique  models  the  effect  of  collisions  which  would 

prevent  the  slower  particles  from  running  away  in  a real  system.  In  the 

runs  investigated  we  have  considered  two  different  values  of  the  beam  to 

background  density  ratio,  nb/nQ  = 1/9  and  1/180.  The  system  length  is 

256  Xpg,  where  XDq  is  the  Debye  length  at  t • 0;  256  beam  electrons  and 

2304  bulk  electrons  are  followed.  Appropriate  measures  are  taken  to 

handle  the  discreteness  of  the  beam  electrons.  Typical  parameter  values 

are  as  follows.  The  beam  momentum  at  t = 0 is  pfe  ■ 7.63  ovTQ,  where  m is 

the  electron  mass.  The  limiting  value  of  the  momentum  above  which  is 

applied  is  6.0  ®vTq-  The  speed  of  light  c is  taken  as  9.0  vTq,  and 
* ®^DC 

Eqc  * moTv — = 0-1»  w^ere  up  the  electron  plasma  frequency  and  e the 
electron  charge.  Each  electron  is  represented  by  a Gaussian- shaped  finite 
size  particle  of  width  1.0  X^q.  A typical  run  extends  to  t * 103w  _1. 

The  early  evolution  of  the  system  is  dominated  by  the  two-stream 
instability  triggered  by  the  low  density  runaway  beam  as  shown  in  Fig.  1(a) 
During  the  early  phase  one  observes  the  exponential  growth  characteristic 
of  such  an  instability.  Unlike  the  E^a  0 case,  the  instability  does  not 
saturate  by  particle  trapping10,  but  rather  the  wave  energy  grows  secularly 
eventually  reaching  a value  roughly  10  times  larger  than  the  non-runaway 
beam  saturation  level,  as  illustrated  in  Fig. 1(a) .This  figure  exhibits  the 
time  evolution  of  the  total  wave  energy  for  n^/ng  » 1/180;  the  continuous 


* 0. 


curve  corresponds  to  ED(-,  = 0.1  and  the  smaller  dotted  curve  to  E| 

The  latter  curve  is  the  well-known  trapped  particle  saturation 
pattern.  The  reason  for  the  enhanced  excitation  of  plasma  waves  is  that 
the  free  energy  of  the  beam  continues  to  be  replenished  through  the  work 
done  by  the  DC  field,  hence  the  beam  nonlinearities  act  as  a converter  of 
DC  to  AC  energy.  Together  with  the  secular  increase  in  the  total  wave  energy 
one  observes  that  the  wave  spectrum  cascades  toward  smaller  wavenumbers  k 
eventually  reaching  k < u^/c. 

An  examination  of  phase  space  (not  shown)  indicates  that  individual 
runaway  electrons  experience  large  momentum  oscillations  due  to  the  inter- 
action with  the  large  amplitude  waves.  In  this  manner  the  spatially 
averaged  momentum  distribution  f(p)  develops  an  elongated  tail  which 
proceeds  to  expand  in  time  toward  large  momenta.  An  example  of  such  a 
distribution  is  shown  in  Fig.  1(b)  for  the  case  n^/^  = 1/9  and  at  <*>pt  = 550. 
As  seen  in  Fig.  1(b),  the  elongated  tail  is  not  smooth,  but  rather  it 
exhibits  the  transient  formation  and  subsequent  destruction  of  small  bumps 
(beams).11  These  bumps, which  act  as  converters  of  DC  energy  into  wave 
energy,  produce  the  rapid  growth  and  damping  of  the  wave  energy  seen  in 
Fig.  1(a)  over  the  interval  Wpt  > 500.  The  measured  relaxation  oscilla- 
tion period  is  56  Wp-1  for  n^/nQ  = 1/180.  This  relaxation  period  may  be 
evaluated  as  follows.  The  Landau  growth  rate  is  given  by 


yL 


n/2 


wj  3f  (p) 

“5p~ 


vs<uk/k  . 


(1) 


Here  3f/3  may  be  estimated  as  nR/  (p^Pj^) » where  Pj^  ( the  maximum 
momentum  cut  off,  p^,  the  effective  gradient  length  of  the  bump,  and 
nR/nQ  are  90  nrvT0,  3mvT0  and  0.016  respectively  from  the  simulations. 
Here  nR  is  the  number  of  runaway  electrons. 
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The  summation  is  taken  over  bumps  in  resonance  with  the  dominant  components 

of  the  wave  spectrum  (modes  4 and  5),  which  yields  y.  - 1 * 10”2u>  . A much 

l p 

smaller  value  would  be  obtained  from  the  overall  (negative)  slope  of  the 
runaway  tail.  Since  the  wave  energy  varies  by  40%  over  the  trough  value, 
this  leads  to  a relaxation  oscillation  period  - 2 — = 54  “p"1*  In 
the  case  n^/n0  ■ 1/9,  the  measured  relaxation  oscillation  period  is  24  “p"1# 
while  an  argument  similar  to  the  above  yields  xL  - 21  u>p-1.  The  relaxation 
oscillation  period  then  increases  as  the  beam  to  plasma  density  ratio  is 
lowered,  as  expected  from  the  growth  rate  expression,  Eq.(l).  The  above 
estimates  also  imply  that  an  argument  based  on  a steady  state  picture  might 
be  inadequate. 

An  assumption  which  permeates  the  quasilinear  description  of  runaway 
behavior  is  that  the  real  part  of  the  dispersion  relation  is  determined  by 
the  background  (slow)  particles.  However,  as  seen  in  Fig.  i(b),the  self- 
consistent  distribution  consists  of  an  elongated  tail.  The  dispersion  relation 
for  such  a distribution  may  be  written  in  the  highly  relativistic  limit  as 

<*)  2(i-n)  nu  2 

1 " <*$-3k2v2Q  p^  w(u-kc)  ' °* 

where  n = nR/(nQ+nR)  and  pM  is  the  maximum  (cut-off)  momentum,  pM  » me. 

The  modes  thus  introduced  may  be  called  an  electroacoustic  branch  due  to 
the  hot  electron  population  (elongated  tail) . The  appearance  of  a disper- 
sionless low  frequency  branch  is  indicative  of  the  Van  Kampen  effect,  since 
the  cut-offs  in  runaway  electron  momenta  (at  0 and  p^)  introduce  the  non- 
cancellation of  the  phases  of  the  Van  Kampen  modes.12  In  Fig.  2 we  compare 
the  prediction  of  Eq.(2)  with  p^/mc  ■ 13.3  and  n ■ 0.016  (solid  curve) 
with  the  simulation  results  (solid  dots)  obtained  by  standard  correlation 
methods.  It  is  clearly13  seen  that  the 


r — — 

a lf.-  ~ ' ' — 'H| 

I 
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elongated  tail  strongly  modifies  the  wave  dispersion,  as  expected.  Further- 
more, the  strong  nonlinearity  produced  by  the  large  amplitude  waves  gives 
rise  to  multiple-harmonic  branches  (crosses)?4  A salient  feature  in  Fig.  2 
is  that  the  system  supports  modes  whose  phase  velocity  is  slower  than  the 
edge  velocity  of  the  expanding  distribution  for  k < oj^/c.  Consequently,  the 
runaway  particles  are  always  in  resonance  with  a collective  mode.  The 
excitation  of  such  modes  retards  the  free  fall  expansion  of  the  runaway  tail. 

Although  the  interaction  of  an  individual  runaway  electron  with  the 
large  amplitude  waves  is  quite  complicated,  the  average  behavior  of  the 
runaway  population  is  remarkably  simple,  as  shown  in  Fig.  3.  The  average 
runaway  momentum  p (i.e.,  average  above  the  cut-off  at  6.0  mvTQ)  is  found  to 
increase  linearly  with  time,  p * p + aeE-^t,  where  the  coefficient  a « 2/3 
for  nb/nQ  = 1/9  and  a ~ 1/2  for  n^/^  = 1/180.  This  behavior  can  be  under- 
stood through  a simple  scaling  argument  by  realizing  that  for  n^/n^^  » 1/9 
the  number  of  runaway  electrons  is  increasing  only  slightly  and  that  the  time 
averaged  distribution  function  approximately  has  a self-similar 
triangular  form  given  by 


f (p)  =£_(!.£_)  (0  < p < ip) 

3p  3p 


(3) 


Conservation  of  momentum, (d/dt) (p2/2)  = eEp, with 

'3p 
o 

and  f 3p 

o 

yields 

P = (2/3)  eE^, 

in  excellent  agreement  with  the  simulations  (see  Fig.  3). 


(d/dt)(p2/2)  = (d/dt)  [(1/2)  [3P  f(p)p2dp]  = (3/2)  pp 

* o 

r 3p  _ 

eEp  e e I f(p)  p • E^dp  * eED(,  p 

• n 


(4) 


(5) 


(6) 


(7) 
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In  addition  to  the  classical  runaway  threshold  (i.e.  the  Dreicer 
field15),  we  find  another  threshold  DC  field  beyond  which  the  electrons 
keep  running  away  toward  infinite  momenta  without  any  strong  collective 
phenomena  taking  place.  This  critical  DC  field,  Ecr,  is  determined  from 
the  peak  wave  amplitude16  excited  by  a beam  when  EDC  = 0,  i.e. 

Ecr  = ^irn^mc2  s(l  + s)'5/2]5* 

where  yq  is  the  relativistic  factor  and  s = g2Y0(n^/2n0) 1//3,  yq  * (1  - 02) 

When  Eno  > E , the  waves  fail  to  sustain  oscillations  and  thus  slow  down 
the  particles.  In  the  case  n^/n^  = 1/9 » a strong  two-stream  instability  is 
still  experienced  in  the  simulation  for  EDC  = 1.0,  while  no  two-stream 
instability  nor  any  strong  wave-particle  interactions  occur  for  E^  = 2.0. 

In  fact,  it  is  observed  that  as  ED(,  approaches  Ecr  the  level  of  AC  fluctua- 
tions dramatically  increases  and  then  for  > E it  remains  at  the  thermal 

DC  cr 

level.  Simulation  values  of  y0  * 1.31,  c * 9.0  v^,  n^/no  * 1/9  and  s * 0.21 
(at  t * 0)  in  Eq.(7)  yield  E * 1.25,  confirming  the  simulation  results. 

For  nb/nQ  = 1/180, this  threshold  is  found  to  lie  between  EpC  =0.1  and 
E^,  = 0.5  in  the  simulation,  while  Eq.(7)  predicts  E = 0.3. 

1/L  CX 

The  physical  picture  which  emerges  out  of  this  investigation  is  the 
following.  The  presence  of  a DC  electric  field  EpC  < Ecr  prevents  the  usual 
saturation  of  the  beam-plasma  instability.  Consequently,  large  amplitude 
plasma  waves  are  excited  and  cause  large  momentum  fluctuations  in  the  runaway 
population.  The  strong  nonlocal  mixing  in  momentum  space  produces  an  elongated 
momentum  distribution  which  supports  Van  Kampen  (electroacoustic)  type  of 
modes.  The  self-consistent  excitation  of  these  modes  gives  rise  to  the 
transient  growth  and  decay  of  bumps  in  the  distribution  function  accompanied 
by  relaxation  oscillations  in  the  wave  energy.  The  drag  produced  by  these 
waves  reduces  the  runaway  momentum  increase  below  the  free  fall  value. 
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FIGURE  CAPTIONS 


(a)  Time  evolution  of  the  wave  energy  in  the  presence  of 
runaway  electrons.  With  * 0.1  (solid  curve),  the  AC  wave 
energy  is  sustained  to  a much  higher  level  than  with  E^  « 0 
(dotted  curve).  = 1/180.  (b)  The  time  averaged  momentum 

distribution  function  at  t « 550  shows  the  elongated  tail. 

V"o  ■ 1/9- 

Electroacoustic  dispersion  relation  measured  (solid  dots)  and 
predicted  by  theory  (solid  curve)  for  n^/no  ■ 1/180.  Note  that 
there  appear  second,  third  and  fourth  harmonic  branches  (crosses) 
in  the  dispersion  relation.  The  spectral  intensity  of  the 
harmonic  modes  is  highest  at  the  crossing  of  the  harmonic  and 
electroacoustic  branches. 

Evolution  of  the  average  runaway  momentum.  The  solid  dots 
represent  the  simulation  results,  the  dashed  line 
p = pQ  ♦ 2/3  eEp^t.  The  average  runaway  momentum  first  follows 
the  free  fall  line  and  eventually  asymptotes  to  2/3  of  the  free 


fall  value. 
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